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Abstract 

A new calculation of nuclear enhancement factor, used in estimation of Galac- 
tic diffuse gamma-ray flux from proton-proton interaction in order to take 
account of heavy nuclei included in cosmic-rays and interstellar matter, is 
presented by use of a Monte Carlo simulator, DPMJET-3. A new value of 
1.8-2.0 in the energy range of 6-1000 GeV/nucleon, slightly increasing with 
kinetic energy of projectile cosmic rays, is about 20% larger than previous 
estimates. 
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1. Introduction 

Majority of cosmic gamma-rays observed by EGRET onboard the Comp- 
ton gamma-ray observatory concentrate around the Galactic plane, and they 
are known as Galactic diffuse gamma-rays [H . Gamma-ray point sources near 
the Galactic plane are observed above this diffuse 'background' gamma-rays, 
and the diffuse emission model based on high-energy particle interaction with 
interstellar matter/field is extremely important for the analysis of gamma-ray 
emitting objects in the Universe. 

The dominant component of Galactic diffuse gamma-rays above 100 MeV 
is well explained by decay gamma-rays from neutral pions generated in nu- 
clear interaction of high-energy cosmic-rays and interstellar matter [sj. 
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In calculation of this process, most authors have computed gamma-ray flux 
generated by neutral pion decay which is produced in high-energy proton 
interaction on proton target (see Kamae et al. for the most recent cal- 
culation). Actually, there are heavy nuclei both in cosmic-rays and target 
matter. This fact is taken into account by introducing a nuclear enhance- 
ment factor, eu, to be multiplied to the result of gamma-ray yield assuming 
cosmic-ray protons on proton target only. Previous estimates are as follows: 



~ 2.0 by Stecker [5|, 1.5 by Cavallo and Gould [gI, 1.6 ±0.1 by Stephens and 
Badhwar [2], 1.45 by Dermer Q, 1.52 + max(l, {Tp/ 100 GeV)^-^^) by Mori Q, 



where Tp stands for kinetic energy of cosmic-ray proton. Mori [9j used the 
low-energy value of 1.52 given by Gaisser and Shaefer [l3| derived for calcu- 
lation of cosmic antiproton flux without a factor related only to propagation, 
and took the spectral difference of heavier nuclei into account. 



Later, Simon, Molner and Roesler [12| calculated interstellar secondary 
antiproton flux with a Monte Carlo model, DTUNUC, to study contribu- 
tions to the antiproton flux from p-nucleus, nucleus-p and nucleus-nucleus 
collisions. They gave the energy dependence on the nuclear enhancement 
factor (which they call 'nuclear scaling factor' in their paper) which differs 
signiflcantly in low-energy region {Tp < 10 GeV) from a constant value as 
given by Gaisser and Shaefer [l3]. The situation is similar to the case of 
Galactic diffuse gamma-rays, and we need detailed calculation of nuclear 
effect taking account of energy dependence for the new era of gamma-ray as- 
tronomy with a launch of Fermi (formerly called GLAST) gamma-ray space 
telescope in June 2008 [11]. 

Here we use a Monte Carlo model, DPMJET-3 [l^, which is an up- 
dated version of DTUNUC, to calculate a nuclear enhancement factor for 
the Galactic diffuse gamma-ray 

2. Calculation 

We calculated gamma-ray yield from p-nucleus, nucleus-p and nucleus- 
nucleus collisions using DPMJET-3. It is a Monte Carlo model capable of 
simulating hadron-hadron, hadron-nucleus, nucleus-nucleus, from a few GeV 



^Recently, Huang et al. calculated gamma-ray yield from (p, a) + (interstellar matter) 
interactions using the same Monte Carlo model pj, but their results are given as final 
spectra and the nuclear effect on the gamma-ray production was not explicitly discussed 
in their paper. 
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Table 1: Multiplication factor at T—10 GeV/nucleon. G&S is quoted from ref . [lot . 



Nuclei 


TTlip 

DPMJET-3 


G&S 


DPMJET-3 


G&S 


H {A=l) 


(1) 


(1) 


3.81 


3.57 


He (A=4) 


3.68 


3.57 


14.2 


12.6 


CNO {A=U) 


11.7 


11.4 (N) 


42.5 


40 (N) 


Mg-Si (A=25) 


20.3 


20 (Al) 


73.2 


71 (Al) 


Fe (A=56) 


38.8 


40 


142 


135 



up to the highest cosmic-ray energies, and has been successfully applied to 



the description of hadron production in high-energy collisions [13| . To check 
our calculation with this model, "multiplication factor", rriip and mjn for 



proton and Helium target, respectively, given in Gaisser and Shaefer [10 
was recalculated for T = 10 GeV where T is a kinetic energy per nucleon. 
This is a relative yield of gamma-rays from nucleus-p and nucleus-Helium 
collisions compared with that from p-p collisions, and is defined as rriip = 
{(^ip/ (Tpp){{n^,ip) / {n^,pp)) and similarly for mj^, where (Tip (iTpp) and (n^,ip) 
((^7,pp)) are a production cross section and an average number of gamma- 
rays produced in nucleus-p {p-p) collision, respectively. The results are given 
in Table [H For each entry, 20,000 collisions were generated. Reasonable 
agreement is seen, and in fact, using the flux values at 10 GeV/nucleon, 
nuclear enhancement factor can be calculated as 

= 1 + E "..m + E X = 1-52 (1) 



and is close to 1.52, the value given by Gaisser and Schaefer [lOj, assuming the 
interstellar matter abundance and cosmic-ray composition given by ref. [l^ 
(with the fraction of Helium in interstellar matter r = 0.07). Thus the 
low-energy value of the nuclear enhancement factor changes little by use of 
modern nuclear interaction model. 

Figure [1] shows the energy dependence of multiplication factor. Due to 
limitations of the DPMJET-3 model, calculations below about 5 GeV/nucleon 
produce errors and we stop calculations there. Table [2] shows multiplication 
factors for heavy target nuclei calculated using DPMJET-3. Note the approx- 
imate symmetry between projectile and target nuclei (ex. p-Fe vs. Fe-p). 
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Figure 1: Energy dependence of multplification factor for various combination of projectile- 
target nuclei. 



Table 2: Multiplication factor at T=10 GeV/nucleon for heavy nuclei target. 



Projectile 




Target 




nuclei 


CNO 


MgSi 


Fe 


H (A^l) 


11.6 


20.1 


38.9 


He (A=4) 


42.3 


73.4 


143.1 


CNO {A=U) 


120.6 


204.2 


386.3 


Mg-Si (A=25) 


204.2 


343.1 


628.8 


Fe (A=56) 


384.4 


634.0 


1067 
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Table 3: Primary cosmic-ray spectral parameters |l5|. 



Component 


a 


K 


b 


c 


H {A=l) 


2.74 ±0.01 


14900 ± 600 


2.15 


0.21 


He (A=4) 


2.64 ±0.01 


600 ± 30 


1.25 


0.14 


CNO {A=U) 


2.60 ±0.07 


33.2 ±5 


0.97 


0.01 


Mg-Si (A=25) 


2.79 ±0.08 


34.2 ± 6 


2.14 


0.01 


Fe (A=56) 


2.68 ±0.01 


4.45 ±0.50 


3.07 


0.41 



Next, we explore other elements of nuclear enhancement factor: cosmic- 
ray spectra and composition of interstellar matter. 

The primary cosmic-ray spectra, 0j(T), are summarized in empirical for- 
mulas by Honda, Kajita, Kasahara and Midorikawa (2004) [l5|: 

(j)i{T) = K{T + 6exp(-cyT))-" [m-^s-^sr-^] (2) 

with parameters a, K, b and c tabulated in Table [3], where cosmic-ray nuclei 
are divided into five groups. If we use these spectra, the nuclear enhancement 
factor is calculated to be 1.66 at 10 GeV/nucleon, about 10% larger than that 



by Ref.flO 



As for the composition of interstellar medium, Meyer summarizes elemen- 
tal abundances in local Galactic interstellar media in detail [l^. In his Table 
2 the He/H ratio is0.096|l which is significantly larger than 0.07/0.93 = 0.075 
assumed in Gaisser and Shaefer l^]. This fact leads to a further larger value 



of eM- Numerically, if we use relative abundance of H:He:CNO:NeMgSiS:Fe= 
1 : 0.096 : 1.38 x lO'^ : 2.11 x 10"^ : 3.25 x 10"^ following the compilation by 



Meyer 16|. Then we obtain the nuclear enhance factor at 10 GeV/nucleon: 



eM = 1 ±e(H) ±e(He) ±e(CNO) ±e(NeMgSiS) ±e(Fe) = 1.84 (3) 

where e{i) is the contribution from interstellar nuclei i bombarded by cosmic- 
rays with the composition given by the flux expressed by Eq.(l2]), except p-p 
interaction, calculated in the same manner as in Eq.([l]). This value is about 
20% higher than previous calculations. As discussed above, most of the 



^This ratio is also consistent with the average of abundance in interstellar gas in various 
objects tabulated in Table 21.9 of Ref.[i3], 0.096 ± 0.07 (omitting 30 Dor and Crab which 
shows extra values). 
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Figure 2: Energy dependence of nuclear enhancement factor contributions from each ISM 
component. 



difference resides in abundance of Helium in cosmic-ray flux and interstellar 
medium. Table |4] shows the contribution of each component to the nuclear 
enhancement factor. 

Fig. [2] shows the energy dependence of the nuclear enhancement fac- 
tor including contributions from each ISM abundance. The total factor is 
also tabulated in Table [51 Although it is nearly constant above about 10 
GeV/nucleon, it decrease gradually toward lower energies. However, the cal- 
culation is limited above above 5 GeV/nucleon, again due to the limitation 
of DPMJET-3. 



3. Conclusion 

Nuclear enhancement factor used in models of the Galactic diffuse gamma- 
rays from cosmic-ray and interstellar matter interaction has been calculated 
by use of a Monte Carlo simulator, DPMJET-3. A new value of 1.8-2.0, 
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Table 4: Nuclear enhancement factor decomposed to each component at 10 GeV/nucleon. 









Targ 


et 




Sum 




H 


He 


CNO 


NeMgSiS 


Fe 




Projectile 














H 


1 


0.405 


0.0177 


0.0047 


0.0006 


1.428 


He 


0.203 


0.083 


0.0036 


0.0035 


0.0004 


0.293 


CNO 


0.038 


0.015 


0.0006 


0.0018 


0.0002 


0.055 


MgSi 


0.033 


0.013 


0.0005 


0.0026 


0.0003 


0.049 


Fe 


0.014 


0.006 


0.0002 


0.0021 


0.0002 


0.022 


Sum 


1.288 


0.520 


0.023 


0.0147 


0.0017 


1.845 



Table 5: Energy dependence of nuclear enhancement factor. 



Kinetic energy (GeV/nucleon) 


6.31 7.94 10.0 100 1000 


Nuclear enhancement factor 


1.75 1.84 1.84 1.89 2.00 



slightly increasing with kinetic energy of projectile cosmic rays in the range 
6-1000 GeV/nucleon, is about 20% larger than previous estimates. This 
factor could have a larger energy dependence below 5 GeV/nucleon where 
DPMJET-3 is not applicable (see Simon et al. 12| for the case of antiproton 
production), but we do not have a proper tool to treat that energy region. 
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